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Abstract 

This paper presents a new framework for asynchrony. This has its 
origins in our attempts to better harness the internal decision making 
process of cellular automata (CA). Thus, we show that a max-plus alge¬ 
braic model of asynchrony arises naturally from the CA requirement that 
a cell receives the state of each neighbour before updating. The significant 
result is the existence of a bijective mapping between the asynchronous 
system and the synchronous system classically used to update cellular 
automata. Consequently, although the CA outputs look qualitatively dif¬ 
ferent, when surveyed on “contours” of real time, the asynchronous CA 
replicates the synchronous CA. Moreover, this type of asynchrony is sim¬ 
ple - it is characterised by the underlying network structure of the cells, 
and long-term behaviour is deterministic and periodic due to the linearity 
of max-plus algebra. The findings lead us to proffer max-plus algebra as: 
(i) a more accurate and efficient underlying timing mechanism for models 
of patterns seen in nature, and (ii) a foundation for promising extensions 
and applications. 


1 Introduction 

A cellular automaton (or CA, where we also abbreviate the plural “cellular 
automata” to CA, allowing the context to save confusion) is a discrete dynamical 
system, consisting of an array of identical cells, each possessing a state. The 
states evolve, according to some local rule, in discrete time steps. The first CA 
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models were synchronous, where all cells update once on each time step, and the 
foundations of the study of these CA were laid by Wolfram in the 1980s [l][2]. A 
popular application of such CA is the construction of models of natural growth 
processes such as seashell patterns and snowflakes |3|4]. Figure]!] (right) shows a 
CA pattern typically examined by Wolfram inscribed on a seashell; the similarity 
to the real seashell pattern (on the left) is evident. A natural extension is 



Figure 1: Seashell patterns: left is real, right generated by a CA. Source: |3 


the introduction of asynchronous update times. Indeed, in terms of seashell 
patterns, Gunji demonstrated that different asynchronous update methods yield 
different patterns, leading to their conjecture that asynchrony is intrinsic to 
living systems [5]. 

A preliminary observation of asynchronous cellular automata was made 
in 161, where the authors compared the properties of synchronous CA with 
two types of CA that iterate asynchronously. Subsequent studies employed 
the methods of 6 as special cases to conduct specific studies into asynchronous 
CA [7-11 . Many of these authors attest to asynchrony as being stochastic 
in nature. This is a general viewpoint in light of their applications: such asyn¬ 
chrony relies on continuous time 11 and is also likely to be more robust EMU 


thereby aiding a better description of biological phenomena. For example, given 
a system of coupled cells, the update of cell states depends on a predefined prob¬ 
ability [6,11 . This consequently also led Schonfisch and de Roos to conjecture 


that, while synchronous updating can produce periodic orbits, asynchronous 
systems will only yield patterns that converge to a fixed point or patterns that 
are chaotic 


11 


The argument for asynchronous updating being stochastic has been chal¬ 
lenged by authors such as Cornforth et al in 12 . They claim that such prob¬ 


abilistic updating schemes are used because of the oversimplification of biolog¬ 
ically inspired models. They further argue for mimicking appropriate aspects 
of nature more closely to create better computational models. Thus, the au¬ 
thors have drawn attention to a large class of behaviours of natural processes, 
in which the updating is asynchronous but not stochastic 
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Moreover, underlying a synchronous update scheme is the notion of a ‘global 
clock’, in the sense that it assigns the same update time to all cells. In recent 
years, the disadvantages of synchrony in this context have been exposed 14 


On the other hand, a remedy has been presented from the perspective of parallel 
computing devices, wherein such a distribution of a global signal proves costly. 
Thus, the proposal is to allow cells to determine their own update time through 
local interactions (see 14 and the references therein, including 15 ). Although 
the storage cost is higher, this scheme points towards a more natural form of 
computation. 

In [I], Wolfram explored synchronous CA on a one-dimensional lattice, where 
cells take the Boolean states 1 or 0. The CA state of cell i was dependent on the 
states of three connected cells, called the neighbourhood of i. These three-cell 
neighbourhood CA were termed “elementary cellular automata” (ECA). 

The ECA may be regarded as special cases of random Boolean networks. 
The different types of these networks were first classified by Gershenson in 
2002 [l6b The types considered included asynchronous random Boolean net¬ 
works (ARBNs), in which nodes are selected to be randomly updated at each 
time step, and deterministic ARBNs (DARBNs), where the node to be updated 
is selected deterministically. Gershenson talked of DARBNs as being more ad¬ 
vantageous because of their modelling capabilities, which are more straightfor¬ 
ward than ARBNs that rely on the stochasticity of asynchronous phenomena. 
Gershenson further proposed DARBNs as better representations of the famous 
genetic regulatory networks of Kauffman 17 , as they are asynchronous but do 
not rely on stochastic methods. 

Following on from Gershenson’s idea of using determinism as a more ‘model- 
friendly’ form for asynchrony, a goal of this paper is to exploit this avenue 
by presenting a new, deterministic framework for asynchrony. No matter how 
well it matches the real system, we claim that the essence of many interesting 
and important asynchronous processes is lost by using probabilistic updating 
schemes. 

As a real example, consider Figure [I] again, showing a seashell pattern; it is 
interesting to see that the same pattern (and many other such seashell patterns, 
as well as growth processes such as snowflakes [4]) may be replicated quite 
accurately by a CA model. Traditionally, some difference between the two 
versions - real and CA - would be ascribed to a fault, random or otherwise, such 
that a better approximation may be obtained by adding stochastic asynchrony 
in the CA rule. As alluded to by Cornforth et al in |l2j, such asynchrony 
has tended to simplify these dynamics into a probability (or the like) of cells 
updating their states. Inspired by the ideas of ‘local clocks over global clocks’ 
and 15 , we will get into the heart of the matter and study the pattern 


14 


of exchanges that takes place locally, that is, before any data is transferred 
between cells. 
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1.1 Network Description of the Cellular Automaton Lat¬ 
tice 

Our work views each cell in a one-dimensional lattice as a processor which re¬ 
ceives input from its neighbourhood. Having received this input, the processor 
computes its new state (as a function of the input states), then sends a corre¬ 
sponding output to its connected neighbours. This type of information exchange 
can be represented by the digraph in Figure [2] where each node represents a 
processor and directed arcs between nodes indicate the direction of information 
transfer. 

Figure [2] shows three arcs pointing to each node, indicating that there are 
three processors - therefore three neighbours (including i itself) - sending infor¬ 
mation to each process i. Thus, the neighbourhood size n of each node is n = 3; 
we also refer to such a neighbourhood as an n-neighbourhood (or n-nbhd). The 
figure particularly shows that each processor sends output to itself as well as to 
its left and right neighbours. We refer to this type of network as a regular n- 
nbhd network or simply a regular network if n is understood. Thus, the regular 
3-nbhd network describes the lattice for the ECA, where cells are depicted by 
nodes. For this reason, we use the terms cell and node to mean the same thing. 


. TltXtXtAl tX' 

Figure 2: Digraph representing a regular 3-nbhd network 


It is also assumed that there is a processing delay associated with each nodal 
computation of CA state. Moreover, we incorporate a transmission delay, which 
is the time taken for a CA state to be transmitted to other cells that require 
it. These two parameters are the means by which we obtain asynchrony: a 
divergence from classical (synchronous) CA models since there the computations 
are assumed to occur instantaneously. 


1.2 Contents 


In Section [2j we review the ECA, a classical synchronous model of CA, which 
is later used to present our asynchronous model. In Section [3j we show how 
max-plus algebra provides a natural way to mathematically model asynchronous 
CA. By covering related graph theoretical techniques and known proven results, 
we also show that this asynchronous system is periodic, characterised by the 
connectivity of the underlying network. It will be seen from the depth of theory 
covered that a max-plus algebraic model of asynchrony is more bespoke and 
addresses all the intricacies of the internal dynamics within a cell (not just the 
external). Nevertheless, those already familiar with max-plus algebra may skim 


over Sections 3.1 and |3.2[ noting the few places where we mention the link to our 
system. The “contour plot” is introduced as a framework for this asynchrony in 
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Section [3~3l This is followed by Section [4j which shows the effect of the max-plus 
asynchrony on cellular automata. We finish with concluding remarks. 


2 Synchronous Cellular Automata 

Let Si denote the state of cell i. The index i denotes the position of the cell 
in the one-dimensional Euclidean plane, so that the state of the CA at a given 
time f £ R is represented by the string si(t)s 2 {t) ■ ■ ■ SN{t) or the vector s(t) = 
(si(t), S 2 (t),..., Sjv(f)), where N is the size of the lattice. The ECA assigns a 
symmetrical neighbourhood of three nodes, where node i is contained in its own 
neighbourhood (as in Figure [ 2 ]). Here, the CA rule is a function / : {0, l} 3 —► 
{0,1} given by s,(f + 1) = /(sj_i(f), Si(t), Sj+i(t)). An example of such a CA 
rule is the following. 


j=i +1 

*(*+!)= E Sj(t) mod 2. (1) 

j=i —1 

The rule in equation is named ECA rule 150 by Wolfram [ 2 ] and this is how 
we refer to it throughout the report. 

Consider a regular 3-nbhd network of twenty cells (connected in a ring such 
that end cells are adjacent). If each cell is depicted by a square - coloured if 
Si = 1, clear if s % = 0 - then the output produced after each of 30 iterations 
of equation 0 is shown in Figure [3j For fixed t. the CA state Si(i) • • • S2o(t) 
represents a horizontal line of cell states, and the initial CA state is 

/ n \ / 1 if i = 10 

“ { 0 otherwise. (2) 

Such an output as Figure[3]is referred to as a space-time pattern (or plot). It is 
evident that the absence of transmission and processing delays means that the 
t th update of each cell occurs at the same time as the f th update of every other 
cell, that is, synchronously. 

2.1 State Transition Graph 

Figure [3] displays periodic behaviour. In other words, for each time step t, the 
CA yields a state s(f), which is seen again after a few more time steps. Given 
the initial CA state s(0) and the CA rule /, an orbit of s(0) is the sequence of 
states obtained by applying / on s(0) sequentially. If / is applied k times, we 
represent this as f k ( s(0)) = /(/(• • • /(s(0)))). We define periodic behaviour as 

k times 

follows. 

Definition 2.1. Consider the CA rule / and network size N. Let s (k) = 
f k { s(0)) for all k > 0, where s(k) is a CA state represented by a 1 x IV vector. 
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Figure 3: CA pattern produced by ECA rule 150. Time travels vertically down, 
cells are labelled on the x-axis. 


For some t > 0, if there exists a finite number p £ N such that s (t + p) = s (t), 
then the set of states 


{s(t),s{t + l),...,s{t + p- 1)} 

is called a periodic CA orbit , where p is the CA period of the orbit. 

As an example, consider the underlying network of size 4 as given in Figure[4j 
We give the CA rule in words: the new state of each cell is the sum of the 
states ( mod 2) of its neighbourhood cells on the previous time step, where 
the neighbourhood of cell i comprises those cells whose outgoing arc points to 
i. (Note that this rule is an extension of ECA rule 150 to arbitrary lattices.) 
For small N, as is the case here, it is useful to represent each system state as 
a vertex in a digraph. Thus, there is an arc from CA state s* to CA state s M 
if and only if f(s l ) = s M . The digraph is known as a state transition graph or 
STG for short. For the system in question, the STG is given in Figure [5j Each 
CA state is shown in string form, where the i th digit represents the CA state of 
the * th node. 

To determine the evolution of the CA, we can follow the arcs in Figure [5j It 
can be seen that all initial CA states asymptotically evolve into four periodic 
orbits, represented as circuits in the STG. Two of these (states 0000 and 1001) 
are period-1 orbits (or fixed points in conventional dynamical systems language), 
and two are period-3 orbits; all other states are transient. Such an STG is an 
artefact of a synchronous CA. 
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Figure 4: Size 4 network of CA cells. 
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Figure 5: State transition graph of ECA rule 150 generalised to a lattice of size 
4 (namely the network given in Figure |4j). 


3 Asynchronous Model 

Consider a synchronous cellular automaton. Due to the synchrony, it is possi¬ 
ble to draw horizontal lines in the corresponding space-time pattern such that 
each line represents the update times of all cells at some fixed time. The CA 
is a discrete time dynamical system, so the horizontal lines may be drawn in 
sequence, evenly spaced, as in Figure [6]( a). We call such a space-time plot a 
contour plot , and each horizontal line is referred to as a contour. The contour 
plot may be thought of as a frame on which the CA states are overlaid and 
simulated. Now consider altering these contours so that cells do not necessar¬ 
ily update synchronously. The corresponding contour plot may then look like 
Figure |6](b), which shows the contours having variable shapes - updates occur 
asynchronously. We shall return to the contour plot after looking at how our 
asynchronous model may produce it. 

We present the asynchronous model as follows. Consider a cell i contained 
in a regular n-nbhd network of N cells. The cell carries a CA state (1 or 0) 
which changes with time depending on the rules that we employ. Thus, we can 
plot points on the real line corresponding to when these changes occur. The real 


7 


































(a) Synchronous contour plot 


(b) Asynchronous contour plot 




Figure 6: (a) Synchronous and (b) asynchronous contour plots. The contours 
indicate update times of cells and act as a frame on which the CA may be 
evolved. Both lattices are connected as a regular 3-nbhd network of 20 cells. 
Time is on the vertical axis and the horizontal axis represents the cell positions. 
In (a), contours are horizontal. In (b), update k of all cells is represented by 
contour k (counting from the top). 


line represents time and the points are the update times of the CA state. Let 
Xiik) denote the k th update time for cell i. We also refer to k as a cycle number. 
Once each cell in the neighbourhood of cell i has completed its fc th cycle, it 
sends the updated state to i. The transmission of such a state from cell j to 
i takes transmission time Tij(k). The update of cell i takes a processing time 
and it is represented in the k th cycle by £,(&). If n = 3, we have the following 
iterative scheme for the (k + l) th update time of cell i. 

Xi(k + 1) = max{ Xi-i{k) + n t i-\{k),Xi{k) + r M (fc), 

x i+1 (k) + T iti+ 1 (k)} + £i(k + 1) (3) 

The above sequence of interactions yielding a state change is depicted in Fig¬ 
ure [7J Notice that we have now expanded on the simpler notion of inter-cellular 
communication by focussing closer on intra-cellular communication, that is, the 
communication within a cell itself. This is the key to our study of asynchrony, 
and it has arisen naturally from the requirement of a CA cell knowing the states 
of its neighbours. 

We refer to these internal processes as events. In Figure [TJ there are two 
significant types of events: “receive” and “send”. The three times Xi-i(k), Xi{k) 
and Xi + \{k) are “send” event times, (i.e., when the corresponding CA states are 
sent). The time 

max{xj_i(fc) + Tj i j_i(fc), Xi(k) + T iti (k), x i+1 (k) + T iti+ \(k)} 

is when node i receives the aforementioned “send” states; it is therefore a “re¬ 
ceive” event. Once received, node i processes its new CA state (by applying 
a CA rule on the received states); this takes time duration ^(fc + 1). Once 
processed, node i sends its state to connected nodes at time Xi(k + 1); this is 
another “send” event. Note that we make no distinction between update times 


































time x,_ x {k) 


time x t {k) 


!,(*+!) 


time x l+l (k) 



T,._i ,(£ + 1) 

A 


send state update time x,.(fc + l ) 


x,.,.(A' + l) t,- + m (* +1 ) 

T A 


Figure 7: The processes internal to the k th state change at cell i. Real time 
travels vertically down. Arrows indicate the destination of the sent state. Labels 
on arrows indicate the time taken for the corresponding process. 


and events; both are the same. We, however, formally refer to Xi(k + 1) as the 
update time of i, as introduced earlier. 


Note 3.1. This idea of a network depicting events in a space-time diagram is 
not too dissimilar to Wolfram’s causal network in 
delve into this relationship later. 


25 Chapter 9]. We shall 


3.1 Preliminaries: Max-Plus Algebra and Graph Theory 

The max operation enables one to interpret equation ^ in max-plus algebra. 
This is useful because it converts a nonlinear system into a linear system in this 
new algebra, which subsequently shares many important features with conven¬ 
tional linear algebra. Particularly novel applications include the modelling of 
railway network timetables 18 , manufacturing processes 19 , and even cellu¬ 


lar protein production 20 . In presenting max-plus algebra, we borrow most 


notation and terminology from 18 


Define e = — oo and e = 0, and denote by 1 
a, b £ Umax, define operations © and © by 


the set R U{ £ }- For elements 


a © b = max(a, b ) and a © b = a + b. 


The set R max together with the operations © and © is what we refer to as 
max-plus algebra and it is denoted by 7?. max = (Rmax, ©, ©, £, e). £ = — oo is the 
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‘zero’ (i.e., \/x £ M max , £©£=£©£=£ and e®x = x®e = x), whilst e = 0 
is the ‘unit’ element (i.e., \/x £ R max , e®x = x®e = x). 

Rmax is associative and commutative over both operations © and © while © 
is distributive over ©. (In addition, © is idempotent in 7£ max , so that max-plus 
algebra is a commutative and idempotent semiring.) 

In this paper, we remove the dependence on k of the processing and trans¬ 
mission times so that equation (J3| is written 

Xi(k+ 1) = max{.T,;_i(fc) + t m _i, Xi(k) + t m , x i+ \(k) + T i}i+ 1 } + &. (4) 

We can now write this in max-plus algebra: 

Xi(k + 1) = { (t m _i © Xi-i(k )) © (t. i ti © Xi(k)) 

©(7i, i+ i ©a; i+ i(fc))}©& (5) 

We often save space and clarify the presentation by omitting ©, much as in 
conventional algebra. Thus, x®y = xy and equation © can be written as 

Xi(k + 1) = £i{Ti t i-iXi-i(k) © T iti Xi(k) © T»,»+lXi+l(fe)} (6) 

Since © is distributive over ©, we can write this as 

Xi(k + 1 ) = £iT iti -iXi-i(k) © £iT iti Xi(k) © tiTi ti+ ix i+ i(k) ( 7 ) 

To represent a max-plus power, we follow from associativity of © and define, 
for x £ l ma x, 

x® n = f x © x © • • • © x (8) 

'---v--' 

n times 

for all n £ N with n / 0. For n = 0, we define x®° = e (=0). 

Max-plus algebra extends naturally to matrices, and this allows the con¬ 
current modelling of the update times for all nodes. Denote the set of n x m 
matrices with underlying max-plus algebra by R™*™. The sum of matrices 
A, B £ R^, denoted by A © B, is defined by 

[A © B]ij = a,ij © bij (9) 

where a^- = [A].y and = [B]^. In the same vein, for matrices A £ and 
B £ the matrix product A © B is defined by 

i 

[A © B]ij = ^ a ik © b k j = max {a ik + b kj } (10) 

fe =i feefi,...,*} 

For a £ R, the scalar multiple a © A is defined by 

[a © A]ij = a © dij. ( 11 ) 

As in classical matrix manipulation, the max-plus matrix addition © is 
associative and commutative, whilst the matrix product © is associative and 
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distributive with respect to ®; it is usually not commutative. Similarly, the 
operation ® has priority over ®. 

The elements of R™ ax = 1“^ are called vectors. A vector is usually written 
in bold, as in x, whilst the vector with all elements equal to e is called the unit 
vector and is denoted by u. 

We are now able to define matrix-vector products. The product A®x, where 
A € R"*™ and x <E R" is 


[A®x] j = (T)a ifc (g> x k = max {a ik + Xk}. 


( 12 ) 


k =1 


pxn 
^max 5 


Moreover, for the square matrix A £ 
defined by 

A® k = f A <g> A <g> • • • <8> A 


denote the fc th power of A by A® fc , 

(13) 


n times 


for all k € N with k ^ 0. For k = 0, we set A® 0 = f E(n, n), the identity matrix 
whose diagonal elements equal e and all of whose other elements are e. 

Having established the preliminaries above, a system of N such equations as 
([7| can now be given in the form 


x(fc + 1) = P <S) x(fc) 


(14) 


where x(fc) = (xi(k), X 2 {k ),..., XN(k)) T . P is the N x N matrix defined by 
Aj ® T, where 




Ac = 


£ 

6 


\ 


£n / 


and 


T = 


Til 
T21 

£ 

\ T N, 1 


Tl2 £ 
T22 T 23 

£ 

£ 


TlAT 

£ 


TAr-l,AT-2 TJV—1,AT—1 Tjv-I.JV 

£ Tjv,1V-1 Tjvjv / 


A 4 is referred to as the processing matrix and T is the transmission matrix. We 
call equation (14) a max-plus system (of dimension N) where the vector x(fc) 
is the state of the system. P is called the timing dependency matrja0 of the 
network of cells. NB: Here, the term “state” refers to update time and is not 
to be confused with “CA state”. Nevertheless, the context should make this 
distinction clear. 


1 This name (along with the later timing dependency graph) is inspired by its use in another 
novel application of to the timing of digital hardware in 121 . 
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To a network of cells, we associate a digraph , such as in Figure [2] In general, 
we define a digraph as G = ( V , A), consisting of a set V and a set E of ordered 
pairs (a, b ) of V. (Often, we refer to the digraph simply as a graph which, in 
turn, also refers to the “network” of our application.) The elements of V are 
called vertices or nodes and those of E are arcs. 

An arc (a, b ) is also denoted ab, and we refer to an arc aa as a self-loop. For 
the arc ab, a is the start node and b is the end node ; a is also referred to as the 
predecessor node of b whilst b is the successor of a. 

By assigning real numbers (called weights) to the arcs of a graph Q = ( V. , E), 
we obtain a weighted graph. The weighted adjacency matrix of Q over M max is 
the matrix W £ 1^” whose (■ i,j) th entry Wij is non-zero ( 7 ^ e) if and only 
if j is a predecessor of i. We also refer to W as a max-plus adjacency matrix. 
Given W £ R™ a ", we denote the associated network as G(W). Thus, P in 
equation (14) is a max-plus adjacency matrix of the regular 3-nbhd network 
since the network is exactly that shown in Figure [2] augmented with arcweights 
fi,Ti r In fact, since our particular max-plus system concerns update times, we 
call Q(P) the timing dependency graph of the system. The neighbourhood of a 
node may now be defined in terms of the adjacency matrix. 


Definition 3.1. Let W be the max-plus adjacency matrix of a digraph of 
connected cells. The neighbourhood of i is A) = {j\wij 7 ^ e}. 

Definition 3.2. Let p = { 0 , 1 , 02 ,... ,a n } be a sequence of arcs. If there are 
vertices vo, v \,..., v n (not necessarily distinct) such that aj = Vj-±Vj for j = 
1 ,...,n then p is called a walk from vq to v n . A walk for which the aj are 
distinct is called a path. Such a path is said to consist of the nodes vq, V \,..., v n 
and to have length n, which is denoted \p\i = n. 

If v n = Vq, then the path is called a circuit. If the nodes in the circuit are 
all distinct (i.e., Vi 7 ^ Vk for i 7 ^ k), then it is called an elementary circuit. 

We define the weight \p\ w of a path p as the sum of the weights of all arcs 
constituting the path. The average weight of p is For a circuit, we refer to 
this quantity as the average circuit weight. 

Definition 3.3. For a graph Q = (V,E), node j £ V is said to be reachable 
from node i £ V, denoted i j, if there exists a path from i to j. Graph Q is 
strongly connected if i —> j for any two nodes i,j £ V. 


Moreover, matrix A £ is called irreducible if G(A) is strongly con¬ 

nected; if a matrix is not irreducible, it is called reducible. Thus, in our max-plus 
system of equation (141, matrix P is irreducible. 


Definition 3.4. Denote the cyclicity of a graph Q by og. 


• If Q is strongly connected, then erg equals the greatest common divisor of 
the lengths of all elementary circuits in Q. If Q consists of only one node 
without a self-loop, then ag is defined to be one 
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• If Cy is not strongly connected, then erg equals the least common multiple 
of the cyclicities of all maximal strongly connected subgraphs (MSCSs) of 
Q. (See Appendix |A| for definition of MSCSs.) 

Definition 3.5. Let A £ R™ a ” be irreducible. The cyclicity of A, denoted 
cr(A), is defined as the cyclicity of the critical graph of A. 

When the matrix is understood, the cyclicity is also denoted by a. 


3.2 Asymptotic Behaviour of the Max-Plus System 


Let x(0) represent the initial state of all cells. Then we can rewrite equation ( fl4| 
as 


x(/fc + 1) = P <g> P <g> ■ • • P <g> x(0) = P®( fc+1 ) (gi x (0) (15) 


or, equivalently, 

x(fc) = P® fe G x(0). (16) 

Given x(0), the sequence of vectors {x(fc) : k £ N 0 }, obtained by iterating 
equation (16), is referred to as the orbit of x(0). A study of such sequences is 
provided in greater scope in 18 Chapters 3 and 4]. We detail the topics relevant 
for this work, with some well-known results expanded upon in Appendix [B| 


Definition 3.6. Let A £ R” a ”. If A € R max is a scalar and v £ R” ax is a 
vector that contains at least one finite element such that 


A ® v = A G v, 


(17) 


then A is called an eigenvalue of A and v is an eigenvector of A associated with 
eigenvalue A. 


For a system with irreducible P such as equation (14) (or ©) , it turns out 
that there is only one eigenvalue, and it is equal to the maximal average weight 
of elementary circuits in Q (P) (See Appendix [Bj Theorem |B.1| ) . Such circuits 
with maximal average weight are called critical and the critical graph of P is 
the graph consisting only of critical circuits in t/(P). 


Definition 3.7. Let A £ R™ ax . For some k > 0, consider the set of vectors 


x(fc),x(fc + l),x(& + 2),... e R” ax 

where x(fc) = A® k ® x(0) for all k > 0. The set is called a (periodic) regime if 
there exists p, £ R max and a finite number p £ N such that 


x(fc + p) = p ® x(fc). 


The period of the regime is p. 

Suppose the initial vector x(0) is an eigenvector of A. Then x(k + 1) = 
A ig> x(fc) for k > 0, so that the period is one. Thus, larger periods are obtained 
when the system is not initialised to an eigenvector. The remainder of Section [3] 
explores this in more detail. We start with the following crucial theorem of 
max-plus algebra. 
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Theorem 3.1. Let A £ R„iax be an irreducible matrix with eigenvalue A and 
cyclicity er. Then there is a fc* such that 

^(g>(fc+< t) _ g, 


for all k > k+. 

Proof. See [l8j Theorem 3.9]. □ 

Let a = cr(P). For an indication of the asymptotic behaviour of our system, 
we apply this theorem to observe the state at epoch k + er for k > k+: 

x(fc + cr) = P® (fc+CT) ® x(0) 

= A® CT <g> P® k ® x(0) 

= A ® CT g)x(fc) (18) 

where \® a is read as A x er in terms of classical algebra. This guarantees the 
periodic behaviour of the max-plus system, where er is the upper bound on the 
period. More specifically, the period p is dependent on the choice of x(0) and, 
since we have seen that the system must also be periodic with period er, we have 
that p is a factor of er. We saw earlier that if x(0) is not an eigenvector, then 
p > 1; we now also know that such a period p will not be larger than er. 

We can use the above to show that the vectors x(fc) in a regime turn out to 
be eigenvectors of P® a associated with eigenvalue A x er: 


A ®"®x(fc) from ^ bOTe P^+-)®x( 0) 

P®‘ 7 ® P® k ® x(0) 

= P®' T ®x(fc). (19) 

In fact, given x(0) and corresponding period p, vectors in a regime are also 
eigenvectors of P®^; this can be shown in the same way as above. 

We now define a measure for the average delay between consecutive event 
times Xi(k) and Xi{k + 1). 

Definition 3.8. Let {xi(k) : k £ N} be an orbit of Xi(0) in K max - Assuming 
that it exists, the quantity Xi, defined by 


Xi = lim 

K —KX> 


Xi(k) 

k 


is called the cycletime of i. 


For an irreducible system such as ours, the vector \ = (xi,X 2 , ■ ■ ■, Xn) of cy- 
cletimes is unique (See Appendix [B| Theorem B.2). Moreover, the irreducibility 
of P ensures that each element of x is the same - specifically the eigenvalue of 
P, that is, 

Xi(k) 

lim —— 

k—> oo k 


= A 
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for any initial condition x(0) € R". (See Appendix |B| Lemma B.l 


Thus, since x is independent of the initial condition, we relate it to our timing 
dependency graph by calling it the cycletime vector of P. P is irreducible in 
this paper, so we let X = Xi f° r any b and therefore refer to the cycletime vector 
of irreducible P simply as the cycletime of P. 

As a compact summary of this subsection, we have shown that asynchrony 
due to an irreducible max-plus system always leads to periodic behaviour, and it 


is characterised by the circuit(s) in Q (P) with largest average weight. Section 3.3 


addresses the impact of the above theory on the asynchronous time framework 
of the contour plot introduced earlier. 


3.3 The Contour Plot 

Figure [7] can be seen as a Hasse diagram of events. We say that events are 
causally related if they are contained in the same chain. For example, the 
“send” events at times Xi-\(k — 1 ), Xi(k — 1 ) and Xi+\(k — 1 ) are not causally 
related since they are not contained in a chain - there is no path of directed 
arcs connecting any of the three events. Consequently, these three events form 
an antichain. (For a formal definition of a Hasse diagram, including chains and 
antichains, we refer the reader to Appendix[D]) By connecting those elements in 
the same antichain, we obtain a piecewise linear plot of the vector x(fc), which 
we define next. 

Definition 3.9. Consider the vector x(fc). A contour is the plot obtained by 
connecting (i,Xi(k)) to (i + l,Xj+i(A;)) with a straight line for each i, ( i = 

!)•••, N). 

Creating a contour for each k gives a pictorial representation of vectors x(fc) 
as a function of k. We call this a contour plot. Figure [ 6 ] displays the contour 
plots of a size 20 system, where the sequence {x( 0 ), x(l),x( 2 ),...} represents 
the contours (counting k from the top). For this reason, we interchangeably 
refer to the vector x(fc) by “the k th contour” from now on. Between successive 
contours, we can imagine there being drawn the internal processes of those in 
Figure [7J 

Consider an example system with x(0) = u that yields the following periodic 
behaviour: x(fc + 1) = 5 (8) x(fc) for k > 3. The contours for this system would 
represent vectors in the periodic regime {x(fc)|x(fc + l) = 5®x(fc), x(0) = u, k > 
3}. The period of a regime and cyclicity are related by 1 < p < a. Thus, if 
er = 1 in this example, then p = a = 1 , so that no other period can be obtained 
for all initial states x(0). Therefore, each contour in the contour plot has the 
same shape (separated by 5 time units) as k —> oo; we call this a limiting shape 
of the contours or a limiting contour. For a larger period, we obtain a different 
set of limiting contours. In particular, we obtain different limiting contour plots 
for the cases p = 1 , 2 ,..., cr, each dependent on the choice of x(0). We will 
formalise this in Section [4] 

The idea of a limiting shape in contours suggests a change of coordinates: 
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Given the irreducible matrix P E 


zNxN 

^max 


with eigenvalue A, let 


x(fc) = A® fc 0 y(fc). 


( 20 ) 


We can think of A® fe as a diagonal matrix or the product of A® fe and the iden¬ 
tity matrix E(N,N). The advantage of this is that such a diagonal matrix is 
invertible, its inverse being the diagonal matrix wi th diagonal entries equal to 
A® _fc . Using this property, we rearrange equation (20 1 to obtain 


y(fc) = A® k 0 x(k). 


( 21 ) 


In other words, y (k) is the limit to which the vectors x(fc) — A k tend to as 
k —> oo. By studying the asymptotic behaviour of y(fc) itself, we can deduce 
the shape of the limiting contour. 

The original system follows the recurrence relation x(fc + 1) = P <g> x(fc) for 
some x( 0 ) £ 


lj)ax- Substitute equation (20) into this to obtain 
A®( fc+1 ) 0 y(fc- 


1) = p<g>A® fe 0y(fc). 


( 22 ) 


Interpreting A® fe as a diagonal matrix again yields 

y(k + 1) = A®" (fc+1) 0 P 0 A®* 0 y(fc) 
= A ®" 1 (8) P 0 y(fc) 

= P (g> y(fc) 


(23) 


where P = A ® -1 0 P represents the normalised matrix of P, equivalently ob¬ 
tained by subtracting the eigenvalue of P from each of its entries. The commu¬ 
nication graph of P is the same as that for P (but with different arcweights) 
so that P is also irreducible. However, the maximum average circuit weight of 
G{P), hence the eigenvalue of P, is zero. 

Moreover, it can be shown that P and P have the same cyclicity (i.e., a (P) = 
er(P) = a). Theorem 3.1 tells us of the asymptotic behaviour of the powers of 
an irreducible matrix. Apply this to P to obtain P®( fe+IT ) = O ®' 7 0 P® fe = P® fc 
for k > k+. Thus, using equation (23), 


y (k + cr) 


= P®( fc + CT ) 0 y (0) 
= P® fc 0 y(0) 

= y (k). 


(24) 


So the limiting contour y (k) is periodic with period cr(P). In fact, like x(fc), 
y(fc) is periodic with period p, dependent on x(0), where p is a factor of a. Note 
that this period now conforms with the traditional dynamical systems definition 
of a period in that the sequence {y(fc)|fc £ No} is not monotonically increasing, 
in contrast to the original sequence {x(k)\k £ N 0 }. 

Loosely speaking, Section |3.3| has shown that there is no unique shape to 
the limit of a contour plot. This yields an interesting feature of the max- 
plus asynchronous model: the asynchrony is related not only to the timing 
dependency graph but also to the system’s starting point in time. 
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3.4 The Eigenspace in Max-Plus Algebra 

The set of all eigenvectors of A G KJJ,” associated to eigenvalue A is the 
eigenspace of A. For max-plus asynchrony, the importance of the eigenspace 
is demonstrated through its links to the contour plot, as shown next. 

Consider the recurrence relation x(fc + 1) = P (g) x(fc) and its correspond¬ 
ing contour plot. In Section |3j we established that the vectors in a regime 
are eigenvectors of p® CT ( p ). In other words, each contour is an eigenvector of 
P®°’( p ). Taking linear combinations of ei gen vectors enables the construction of 
the eigenspace of P® CT ( p ) (see Appendix O). This eigenspace is the set of all 
possible periodic regimes (i.e., of all periods p, 1 < p < a(P ), obtained for all 
initial states x(0)), which corresponds to the set of all contour plots that can 
be obtained. 

Theorem |C.1| (Appendix [C]) gives a method for constructi ng t he eigenspace 
of an irreducible matrix. Let B = p® <T ( p ). Applying Theorem C.l to irreducible 
B will yield its eigenspace and consequently all possible contour plots for the 
system x(fc + 1) = P ® x(fc). 

Note 3.2. If P is irreducible, then it is not necessarily the case that B = p® <T ( p ) 
will also be irreducible. Nevertheless, in a highly connected lattice such as the 
regular network of this paper, it is more likely that B is irreducible. If B is 
reducible, then other methods to Theorem |C.1| must be employed. One such 


method is Howard’s policy improvement scheme (see [18] or 22 ). 

Lemma 3.1. Consider a max-plus system having irreducible timing dependency 
matrix P with cyclicity a. For p fixed (1 < p < a), period p regimes are not 
necessarily unique. 


Proof. Each contour in the periodic regime is an eigenvector of P® p . By taking 
linear combinations of eigenvectors, it is possible to construct linearly indepen¬ 
dent eigenvectors of P® p such that the corresponding contours are also linearly 
independent. Thus, there is no unique periodic regime of period p. □ 

The top of Figure [6j( a) shows the limiting contours of a size 20 system. The 
contours depict a period 1 regime, therefore have the same shape. Lemma [3.1| 
says that, if a > 1 for that system, then other shapes of limiting contour may 
be possible for p = 1 periodic regimes. In other words, the limiting contour 
plot is not necessarily fixed, despite p.\ being fixed. The significance of this 
is that a corresponding CA is asymptotically not unique, but dependent on 
the initial time x(0); this can affect the time Xi(k) of the k th update at node 
i relative to Xj{k) (at node j), even when the ‘pattern’ of consecutive update 
times (x(fc), x(fc+l),..., x(/c+p), ...) is independent of x(0) (since the cycletime 
is independent of x(0)). 


4 Cellular Automata in Max-Plus Time 

We now present the first formalism for implementing a cellular automaton asyn¬ 
chronously such that update times are determined by a max-plus system. 


17 






4.1 Contour Plot as a Foundation for Cellular Automata 


Let Si(k) denote the CA state of node i at epoch k. We are concerned with 
Boolean CA states, so that Si(k) £ {0,1}. The unit k is as used in the max-plus 
model which updates the times x(fc). Thus, to be precise, Si(k) is the CA state 
of node i at time Xi (k) . The CA state of the system is represented by the string 
s\{k)s 2 {k) • • • s N {k), which can also be read as the state of all nodes on contour 
k. As a consequence, just as we represented the vector x(fc) by a contour, we 
can represent the CA state s(fc) by the same contour but with the addition that 
the coordinates (i,Xi(k)) now display the state Sj(fc) (e.g., in coloured form, 
where two different colours are used to distinguish the two states 1 and 0). 

Recall the main events that are internal to node i; they occurring within 
cycle k. These events are grouped in two: “receive” and “send”. Once the 
“receive” CA states have all arrived, node i applies a CA rule on this set, to 
obtain the new state Si{k). If all nodes have neighbourhood size n, the applied 
CA rule is the function / : {0,1}" —> {0,1} and the new state Si(k) is calculated 
as 

Si(k) = f{sMi{k - !)). (25) 


4.2 Cellular Automaton Space-Time Plot 

The classical one-dimensional CA is synchronous, so that the k th update time of 
each cell is the same. Consequently, we can think of such a system as having a 
contour plot that contains only horizontal contours. Updates of the CA state of 
the system take place every one time unit, thereby giving the synchronous CA 
a cycletime of 1. The time between contours in this system is thus of duration 
one, although no such duration is depicted; for example, if Sj(fc) = 1 for all k, 
then this is shown as a continuous vertical coloured block in position i. 

Despite varying contour shapes dictating the varying time gaps between 
contours, it is simple and intuitive to construct the space-time plot for CA in 
max-plus time. To illustrate, we use the example of a regular 3-nbhd network 
with size N = 10. Let the positive (diagonal) entries in matrix be represented 
by the vector £ of processing times. We choose the entries in f at random with 
equal probability from all integers between 1 and 30, whilst the non-zero entries 
in T are selected likewise from the integers between 1 and 10. Taking the initial 
time x(0) = u, we obtain a contour plot of update times by iterating the max- 
plus system. We now address the CA state by assigning the depicted space 
between contours as memory: for each node, the CA state remains fixed until 
the time of update, which corresponds to a contour. For node i, the time t, that 
elapses between contours implies that the storing of the CA state in memory 
can be represented as a vertical block of length t, (which is coloured accordingly, 
depending on the CA state). Correspondingly, this may be depicted in a space- 
time plot, the construction of which is shown in three stages in Figure [S| 

At this juncture, it is important to distinguish between the variables t and k. 
The term “time” (or “real time”) now refers to a point t £ M + ; it can be thought 
of as time as we know it. k maintains its role as a discrete epoch. Thus, node i 
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carries a CA state for every point in real time. If we denote the state of node i 
at real time t as s , the CA state of node i can now be understood in two ways: 
Si(k) denotes the CA state on contour k, and it is discretely dynamic, whilst sf ' 1 
represents the state in a dynamical system with a continuous underlying real 
time t. Thus, for example, on contour k = 2, if t = 5.8, we have Si(2) = s[ 5 ' S \ 




(c) (d) 


Figure 8: Construction of the CA space-time plot in max-plus time. The CA 
rule is ECA rule 150. The initial CA state on contour 0 is S5(0) = 1, Si(0) = 0 
for 5. State 0 is coloured dark, state 1 is light. In all figures, the vertical axis 
denotes real time, travelling down, (a) Contour plot with CA states indicated 
on each contour, (b) CA states indicated for all time by filling spaces between 
contours with memory, (c) Contours and space between nodes removed to obtain 
the CA space-time plot, (d) Classical (synchronous) CA space-time plot. 

Figure [8](d) is for comparison with Figure |s](c) and it shows the classical 
synchronous CA having the same initial time x(0), initial CA state s(0) and 
CA rule. Whilst sharing initial conditions and CA rule, the state sW in both 
patterns will generally differ. (This can be seen by simply drawing a horizontal 
line across both patterns at time t and reading off the state of each node at 
that time). The difference in pattern is clearly due to the asynchrony of the 
max-plus system in Figure [8](c) . We show next that we can characterise this 
difference somewhat and, in fact, map the synchronous CA to the max-plus CA 
via the contour plot. 
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4.3 Bijection 

Recall the state transition graph. In Figure |8jd), the CA period is 6 and the 
periodic orbit is the following set. 

{ 0000100000 , 0001110000 , 0010101000 , 0110101100 , . . 

1000100010 , 1101110110 }. 

If we consider the CA states on only the contours in Figure |8jc) (which is seen 
better in Figure (8(a)), we see that they are exactly the same as the states in 
Figure (8]]d). This is a consequence of the max-plus model requiring all neigh¬ 
bourhood states to arrive before processing new CA states. Indeed, this notion 
has also been mentioned (albeit briefly) in (25, Page 1035] under the heading of 
“Intrinsic synchronization in cellular automata”. 

Given the same initial CA state s(0) and CA rule, let S and M. denote 
the orbit of s(0) generated in the synchronous system and the max-plus system 
respectively. Let s g(k) denote the CA state after k iterations of the synchronous 
system; s^v-t (k) denotes the CA state after k iterations of the max-plus system. 
The model uses the same CA rule, applied to the same neighbourhoods, the only 
difference being that the time of application of the rule is different. Then, after 
k iterations of both systems, we clearly have s$(k) = Sj^(k). This defines a 
one-to-one and onto mapping - a bijection - between S and M 1 and we say that 
both systems have the same state transition graph (defined as the transitions 
between states on contours). Thus, the max-plus system need not evolve the 
CA concurrently since the CA plot for the max-plus system may be obtained 
from this mapping. 

In summary, the STG provides a deterministic form for predicting the be¬ 
haviour of the CA in max-plus time. Each state SM(k) in the STG does not 
necessarily correspond to a state in real time t £ R. due to contours not neces¬ 
sarily being horizontal. Nevertheless, the property of memory can be applied to 
ascertain such real time states. 


5 Conclusion 

We have shown that, when modelling discrete asynchronous systems, more at¬ 
tention needs to be paid to the internal processes of a cell. This has resulted in 
the uncovering of a useful theory - that of max-plus algebra. 

Cellular automata are naturally well-suited to be modelled in max-plus time 
because this model requires update on knowledge of all neighbours. This has 
the additional benefit of cells updating only when they are ready. Thus, whereas 
a fixed, periodic, global update time (as in the synchronous case) can be slower 
and less energy efficient, a max-plus asynchronous model consumes only the 
time and energy that local neighbourhoods require. 

The classical ECA corresponds to a strongly connected network in our model. 
Thus, P is irreducible, and this ensures periodic behaviour, which is not what 
is usually associated with the word “asynchrony”. We comment that periodic 
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behaviour is expected also when P is reducible; in this case, the theory is similar 
to what we have covered here but the resultant cycletime vector is not necessarily 
uniform 18 . We have therefore seen that not only is it an efficient system for the 


timing of asynchronous CA, the max-plus system is also a simple, deterministic 
asynchronous model. 

We can get a visual sense of such periodicity via the contour plot, which 
couples the asynchronous update times with the CA. The eigenspace of P® <T ( p ) 
yields the range of contour plots that can be generated due to P; this enables 
us to get a sense of the range of the corresponding CA space-time patterns 
that are possible. By identifying a bijection between synchronous and max-plus 
CA, we can further narrow the aforementioned range of CA patterns since CA 
states on contours can be obtained directly from the synchronous CA. We note 
that this notion of memory has been suggested previously (see [hi], 15 , 25 


and the references therein); this bijection is captured under the topic of “causal 
invariance” in ]25] , wherein the the same causal network (i.e., structure of events 
- past and present) emerges, irrespective of CA states. In this paper, however, 
we have identified a numerical link between the idea and the theory in the form 
of max-plus algebra. 

Importantly, the contour plot also allows one to see that the system visits 
many more ‘interim’ states (in memory) in real time - these are transient CA 
states that illustrate the local dynamics. Thus, although the STG is the same 
as the synchronous STG, max-plus CA gives more information: we can now 
visualise exactly how one state in the STG evolves into another. Max-plus 
algebra offers promising scope for assigning numerical measures to the states 
in memory; the parameters G, 7»j, along with the cycletime, are likely to play 
major roles for this purpose. Such classification scales up to the CA space-time 
pattern as a whole, and attempts to do this are the focus of further work. 

Synchronous CA have previously been held up as models for patterns seen in 
nature (such as those seen on seashells and the growth of snowflakes). Certain 
probabilistic CA improved on this to account for random fluctuations in the 
growth processes. We hypothesise that a max-plus algebraic approach adds 
further realism because it considers the actual processing and delay times that 
may be present within the chemical reactions and relies not on probability. One 
might conjecture the aforementioned transient states in memory to be precisely 
the processes observed during the natural construction of such patterns as found 
on seashells and snowflakes, that is, apparent deviations from ‘normal’ growth 
may, in fact, be part of a transient phase as opposed to some random fault. 
Noting the absence of significant literature on the matter, we believe Boolean 
networks provide a notable avenue to exploit; it would be interesting to see how 
the various topologies of these networks behave under max-plus time. How does 
a max-plus model impact on applications such as Kauffman’s genetic regulatory 
network in 17 ? Moreover, in light of connections to the idea of causal invariance 
in [25| , might max-plus algebra now provide a numerical measure for the related 
topics (e.g., of space, time and relativity) contained therein? 

In addition to the above work, this paper has also laid the theoretical ground¬ 
work for extending the asynchronous model itself. For instance, one such exten- 
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sion includes the minimum operator, and can subsequently be used to conduct 
further studies of asynchronous CA, particularly to better describe an intended 
application. Preliminary results of these new “max-min-plus” models are less 
predictable, particularly manifested by an absence of the bijection described in 


Section 4.3 23 
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A Maximal strongly connected subgraph 

We say that node i communicates with node j, denoted i £A j, if either i = j 
or i —¥ j and j —> i. Note that we allow a solitary node to communicate with 
itself, even if there is no self-loop attached to it. 

It is, thus, possible to partition the node set V of a graph into disjoint 
subsets Vi such that V = V\ U V2 U- • • U V q , where each subset Vi contains nodes 
that communicate with each other but not with other nodes of V. By taking V t 
together with arc set Ei, each of whose arcs has start node and end node in V/, 
we obtain the subgraph Qi = (Vi, Ei). We call this subgraph a maximal strongly 
connected subgraph (MSCS) of Q = (V,E). 


B Asymptotic behaviour of the max-plus sys¬ 
tem 


Theorem B.l. Let A £ be irreducible. Then A possesses a unique 

eigenvalue, denoted A(A), which is finite e). Moreover, this eigenvalue is 
equal to the maximal average weight of elementary circuits in Q(A). Let c 
denote an elementary circuit of Q (^4). Denote the set of all elementary circuits 
of g(A) by C(A). Then 

K A ) = maxjp. (27) 

ceC(A) |c|z 


Proof. See 18 Theorem 2.9]. 


□ 


The following theorem shows that the cycletime vector 


X' = 



xi(k) 

k 


lim 

k —^00 


x 2 (k) 

k 


lim 

k—yoo 


x N (k) 

k 


T 


is unique. 

Theorem B.2. Consider the recurrence relation x(/c + 1) = A®x(k) for k > 0 
and A £ R”);" irreducible. For some x*(0) £ R„ lax whose elements are all finite, 
if the limit lim^oo A exists, then this limit is the same for any initial 

condition x(0) £ R^ ax whose elements are all finite. 
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Proof. See 18 Theorem 3.11]. 


D 


In fact, the condition of irreducibility can be relaxed to that of reducibility, 
as long as each node has at least one predecessor node in G(A); this corresponds 
to all rows of A containing at least one non-zero — oo) element. Thus, in 
the latter (relaxed) case, we obtain a cycletime vector whose elements may not 
necessarily be identical. However, for A irreducible, each element of \ turns 
out to be the same - specifically the eigenvalue of A - as stated in the following 
lemma. 


Lemma B.l. For the recurrence relation x(fc + 1) = A ® x(fc) with k > 0, 
let A £ ®ma" be an irreducible matrix having eigenvalue lei. Then, for 


i = 1,2,... ,n, 


lim 

k—too 


Xi(k) 

k 


= A 


for any initial condition x(0) £ R”. 


Proof. See 18 Lemma 3.12]. 
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C Eigenspace of an irreducible matrix 


As in conventional linear algebra, eigenvectors are not unique in max-plus al¬ 
gebra because they are defined up to scalar multiplication. (It can easily be 
shown that, if v and w are eigenvectors of A £ associated with eigenvalue 

A, then, for a,f3 £ R max , a®v®/3®wis also an eigenvector of A.) 

Consider the definition of the Kleene star for any A £ R™*": 

OO 

T"0#'. (28) 

fc =0 

It is known that, if circuit weights in G(A) are nonpositive, then the Kleene 
star of a square matrix over R max exists 181. Denote the critical graph of A as 
G cr (A) = (V cr (A), E cr {A )) and the normalised matrix A = —A ® A. 

Theorem C.l. Let A £ R]^*™ be irreducible and consider A* to be the Kleene 
star of A = —A ® A. 

1. If node i belongs to Q cr (A), then [A*].i is an eigenvector of A. 

2. The eigenspace of A is 

V{A) = < veR” ax |v= (0 aj® [A*].j for a* £ R max > . (29) 

I i&V cr (A) I 


3. For i,j belonging to G cr (A ), there exists a € R such that 

a® [A*].i = [A*].,- 

if and only if i and j belong to the same MSCS of G cr (A). 


Proof. In 18 Thoerem 4.5]. 


(30) 

□ 
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D Hasse diagram 

With elements taken from [24], we define the Hasse diagram formally. We first 
require the “happened before” relation, denoted by 

Definition D.l. The relation -< on a set of events is defined by the following 
conditions: 

1. If the events a and b are processed by the same processor, and a occurs 
before b, then a b. 

2. If a is the sending of a message by processor A and b is the receipt of the 
message by another processor B , then a -< b. 

3. If a -< b, and b c, then a -< c. 

We say that two distinct events a and b are concurrent if a b and b A a. 
We also assume the properties of irreflexivity, that is, s/s, and antisymmetry 
on the times of events, that is, if a* and bt represent the times of events a and 
b, then a t b t , b t A a t =$_a t = b t . 

In fact, Definition |D.1[ along with the properties of irreflexivity and anti¬ 
symmetry, defines a partial ordering on the set X of all events in our system. 
We say that the set X, along with the relation -<, forms a partially ordered set 
(“poset”). In the following, the pair (X, -<) denotes a poset. 

Definition D.2. Let x and y be distinct elements of a poset (X, -<). y is said 
to cover x if x -< y but no element z satisfies x ^ z < y. 

Definition D.3. Let the set X-\ of n elements {xi, X 2 , ■ ■ ■ , x n } be a subset 
of (X, -<) such that each element may be totally ordered according to -< as 
Xi -< X 2 -< ■ ■ ■ -< x n . Then X\ is a chain. The subset X 2 £ A' is called an 
antichain if and only if no elements of X 2 may be totally or partially ordered. 

These definitions now allow us to formulate the following definition. 

Definition D.4 (Hasse diagram). The Hasse diagram of a poset (A, -<) is a 
graph drawn in the Euclidean plane such that each element of the poset is 
represented by a unique vertex in the graph. Each covering pair x -< y, is 
depicted by a directed arc from x to j/, where the point representing x is below 
the point representing y (i.e., it has smaller T-coordinate). 
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